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ABSTRACT 

This paper explores the phenomenon of energy relaxation for stars in a galaxy em- 
bedded in a high density environment that is subjected continually to perturbations 
reflecting the presence of other nearby galaxies and/or incoherent internal pulsations. 
The analysis is similar to earlier analyses of energy relaxation induced by binary en- 
counters between nearby stars and between stars and giant molecular clouds in that 
the perturbations are idealised as a sum of near-random events which can be modeled 
as diffusion and dynamical friction. However, the analysis differs in one important 
respect: because the time scale associated with these perturbations need not be short 
compared with the characteristic dynamical time to for stars in the original galaxy, 
the diffusion process cannot be modeled as resulting from a sequence of instantaneous 
kicks, i.e., white noise. Instead, the diffusion is modeled as resulting from random kicks 
of finite duration, i.e., coloured noise characterised by a nonzero autocorrelation time 
t c . A detailed analysis of coloured noise generated by sampling an Ornstein-Uhlenbeck 
process leads to a simpling scaling in terms of t c and an effective diffusion constant 
D. Interpreting D and t c following early work by Chandrasekhar (1941) (the 'near- 
est neighbour approximation') implies that, for realistic choices of parameter values, 
energy relaxation associated with an external environment and/or internal pulsations 
could be important on times short compared with the age of the Universe. 

Key words: galaxies: structure - galaxies: kinematics and dynamics 



1 MOTIVATION 

Dating back to the pioneering work of Chandrasekhar in the 
1940's (cf. Chandrasekhar 1943a), it has been recognised 
that discreteness effects will induce changes in the energy 
(and any other collisionless invariants) in a system of stars 
idealised in a first approximation as a collisionless equilib- 
rium. In particular, by modeling interactions between indi- 
vidual stars as an incoherent sum of binary encounters (grav- 
itational Rutherford scattering), Chandrasekhar computed 
a characteristic relaxation time ta on which the energy of 
a typical star will change by a factor of order unity. For 
small systems like open clusters and globular clusters, tn is 
typically short compared with tn, the age of the Universe, 
and this energy relaxation has significant observable conse- 
quences, e.g., for the development of core-halo structures in 
globular clusters {cf. Binney & Tremaine 1987). 

By contrast, for systems as large as entire galaxies, tn is 
typically long compared with tn, so that it has become cus- 
tomary (c/. Binney & Tremaine 1987) to ignore interactions 
between individual stars. However, such interactions are not 
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the only irregularities which can induce energy relaxation. 
For example, in the 1950's Spitzer and Schwarzschild (1951, 
1953) argued that interactions of stars with gas clouds in 
the Galactic disc will tend to increase their random veloci- 
ties and their distance from the galactic plane on time scales 
short compared with tn, thus providing an explanation of 
the observed correlations between kinematics and spectral 
type for Population I stars {cf. Mihalas & Binney 1981). 

The basic physical idea underlying this earlier work is 
that interactions between stars and/or gas clouds can be 
idealised as a sum of instantaneous near-random events, re- 
sulting in a phase space diffusion where, at least in a first 
approximation, the root mean squared change in quanti- 
ties like the energy grows as the square root of time, i.e., 
SE rma oc t 1 / 2 . In particular, in his classic review article on 
'Stochastic Problems in Physics and Astronomy,' Chan- 
drasekhar (1943b) recast the problem in the language of 
nonequilibrium statistical mechanics, allowing for the com- 
bined effects of diffusion and dynamical friction, with the 
effects of diffusion being modeled as Gaussian white noise, 
i.e., a sequence of instantaneous random impulses. 

The objective of this paper is to argue that there are at 
least two other sources of near-random irregularities which, 
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under certain circumstances, can induce significant energy 
relaxation on time scales short compared with tn- 

(1) Real galaxies are all perturbed to a lesser or greater 
extent by other nearby galaxies and, under appropriate cir- 
cumstances, one might expect that these perturbations can 
be approximated as a superposition of near-random inter- 
actions of finite duration, i.e., coloured noise. If the galaxy 
be situated in a comparatively low density environment, the 
dominant sources of irregularities may be a small number of 
satellite galaxies or, perhaps, a single large neighbour which 
execute near-periodic motions relative to the original galaxy. 
In this case, it is clearly inappropriate to approximate the 
perturbations as random events. However, a galaxy embed- 
ded in a high density rich cluster will be exposed to a more 
irregular set of perturbations, associated over the course of 
time with a larger number of different galaxies, and it would 
not seem unreasonable to suppose that these perturbations 
could be modeled as an incoherent sum of near-random in- 
teractions. 

(2) There is also growing evidence, both numerical and ob- 
servational, that real galaxies can exhibit finite amplitude 
oscillations which persist for times long compared with a 
characteristic dynamical time, to- If these oscillations are 
dominated by a small number of normal or pseudo-normal 
modes, they should be well approximated as nearly peri- 
odic perturbations. If, however, the oscillations involve a 
comparatively large number of different modes, it would 
seem that their effects would be better modeled as a sum 
of near-random perturbations, involving a superposition of 
a large number of different frequencies. A collection of dif- 
ferent forces characterised by a random combination of fre- 
quencies combined with random phases is equivalent math- 
ematically to (in general) coloured noise (cf. van Kampen 
1981), with a finite autocorrelation time. 

Section 2 begins by formulating the problem of energy 
relaxation abstractly in terms of D, the diffusion constant, 
r), the coefficient of dynamical friction, and t c , the auto- 
correlation time, i.e., the characteristic time scale on which 
the random irregularities change. These quantities are then 
related to (1) the physical properties of a rich cluster, by 
assuming that the random perturbations acting on a given 
galaxy at given instant are associated with the effects of 
one or two particularly proximate galaxies (Chandrasekhar's 
[1941] 'nearest neighbour approximation'); and (2) the ef- 
fects of incoherent internal pulsations, by exploiting simple 
dimensional arguments. Section 3 summarises a numerical 
investigation of energy relaxation in which diffusion is mod- 
eled as coloured noise sampling an Ornstein-Uhlenbeck pro- 
cess, demonstrating the existence of a simple scaling in terms 
of D and t c . Section 4 concludes by showing that, for reason- 
able choices of parameter values, energy relaxation induced 
by external perturbations and/or incoherent internal oscil- 
lations can be significant over the age of the Universe, and 
then speculating on possible implications for real galaxies. 



2 THEORETICAL CONSIDERATIONS 
2.1 The abstract problem 

Begin by formulating the problem of energy relaxation ab- 
stractly for an ensemble of point masses moving in a time- 
independent potential V(r), each of which is also subjected 



to dynamical friction and coloured noise, i.e., random kicks 
of finite duration. This leads naturally to the consideration 
of a Langevin equation of the form (cf. van Kampen 1981) 



d 2 x a 



<9U(x) 



(a = x,y,z) 



(1) 



dt 2 dx a 

where rj is the coefficient of dynamical friction and F is the 
random force, which is treated as a stochastic variable. As- 
suming in the usual fashion that F corresponds to homoge- 
neous Gaussian noise, its statistical properties are charac- 
terised completely by its first two moments, namely 

{F a {t)) = and 

{F a (ti)F b (t 2 ))=5 ab K(ti-t 2 ), (a,b = x,y,z). (2) 

For the special case of white noise, seemingly appropri- 
ate, e.g., for modeling gravitational Rutherford scattering of 
nearby stars, the autocorrelation function K is proportional 
to a Dirac delta, so that 



K(r) = 2 V QS D (r) 



(3) 



Here the normalisation ensures that the friction and noise 
are related by a Fluctuation-Dissipation Theorem in terms 
of a 'temperature' 0. It is well known that such a Langevin 
description is equivalent (cf. Riskin 1989) to a Fokker-Planck 
equation of the form 



dt cbc 
where 



dV(x) df_ 



<9v 



D 



drK(r) = 267/ 



(4) 



(5) 



is the diffusion constant. 

For coloured noise, the forces at nearby instants of time 
remain correlated so that K is no longer proportional to 
a Dirac delta. As a simple example, one can consider the 
so-called Ornstein-Uhlenbeck process (cf. Uhlenbeck & Orn- 
stein 1930, van Kampen 1981), for which K decreases expo- 
nentially in time. Here 

f/9 



K(r) 



t, 



exp(-M/f c ), 



(6) 



where the autocorrelation time tc sets the time scale on 
which the fluctuating random forces F change appreciably. 
The normalisation entering into this equation ensures that 
the effective diffusion constant D, again defined by eq. (5), 
is fixed completely by r\ and <d. 



2.2 Intuitive expectations 

The effects of energy diffusion induced by friction and white 
noise can be derived directly from the Fokker-Planck de- 
scription associated with the Langevin equation (1). As 
shown, e.g., in Habib, Kandrup, & Mahon (1997), the 
Fokker-Planck equation (4) implies that the mean squared 
change in energy associated with multiple noisy integrations 
of the same initial condition satisfies 



D{v 2 )+ V ((v 2 ) 2 -(v 4 ))+r,((v 2 ){x 2 )-{v 2 x 2 )).(7) 



d((8Ef) 
dt 

However, at early times one can approximate (v 4 ) ~ (v 2 ) 2 
and (v 2 x 2 ) ~ (v 2 )(x 2 ), so that 
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dt 



D(v 2 ). 



(8) 



Noting further that {v 2 } ~ \Eq\, where Eo is the initial en- 
ergy, one then infers that 

d((5E) 2 ) 



dt 



D\E \ 



(9) 



which implies a fractional root mean squared change in en- 
ergy 

Given that an ensemble of different initial conditions each 
evolved in the presence of friction and noise can be viewed 
as a sum of random samplings of Langevin simulations with 
different initial conditions, each satisfying eq. (10), it is clear 
that this relation should also hold quite generally for any 
ensemble of orbits with the same initial energy. The validity 
of this relation has been established numerically for a variety 
of different two- and three-dimensional potentials (cf. Habib, 
Kandrup, & Mahon 1997). 

The obvious question is: how will this result be changed 
if one considers coloured noise with a finite autocorrelation 
time? Recent work indicates that coloured noise impacts or- 
bits via a resonant coupling between the natural frequencies 
associated with the perturbation and the natural frequen- 
cies of the orbits (Pogorelov & Kandrup 1999, Kandrup, 
Pogorelov, & Sideris 2000). White noise is characterised by 
a flat spectral density with power at all frequencies (the 
square of the spectral density is the Fourier transform of the 
autocorrelation function K), so that it can couple to more or 
less anything. Replacing white noise by coloured Ornstein- 
Uhlenbeck noise with a nonzero autocorrelation time yields 
instead a spectral density 



S(u) oc 



(11) 



~> 2 + a 2 ' 

with a — tc , which cuts off for lo 3> ct. Given that the 
characteristic frequencies associated with the unperturbed 
orbits typically scale as t D x , with to a characteristic dy- 
namical time, one might therefore expect (i) that coloured 
noise with t c <^to will have virtually the same effect as 
white noise; but (ii) that coloured noise with the same D 
but t c ^> to will be significantly less efficient than white 
noise in inducing energy relaxation. This in fact turns out 
to be the case. The numerical experiments described in the 
following section lead to the simple scaling 

for t c < to 

(12) 

for t c > to- 



8E r 



\E\ 



Dt 

W\ 



1/2 



1 



2.3 Connection with physical parameters 

It remains to relate D, n, and t c to observable parameters 
characterising a real galaxy, considering first the effects of a 
high density environment. 

Begin by observing that, on dimensional grounds (cf. 
eq. [6]), the diffusion constant 

D ~ F 2 t c , (13) 
where F denotes the characteristic magnitude of the force 



per unit mass associated with the perturbing galaxies. Fol- 
lowing Chandrasekhar (1941) or Chandrasekhar & von Neu- 
mann (1942), suppose now that these effects are associated 
primarily with one or two nearby galaxies. It then follows 
on dimensional grounds that 



Fr- 



m 



and 



Vs 



(14) 



(15) 



where M g represents the mass of a typical perturbing galaxy, 
R s the typical separation between galaxies, and Vs the typ- 
ical relative velocity between galaxies. This implies that 



D 



G Mg 

~mv7 



(16) 



However, by interpreting the 'temperature' O as a charac- 
teristic squared velocity with which the perturbing galaxies 
move through space, one has Q ^ V 2 , so that 



n < 



G Mg 

~mvF 



(17) 



It is easily seen that an application of a similar argu- 
ment to stars moving within a galaxy leads essentially to 
the standard expression for the coefficient of dynamical fric- 
tion associated with binary encounters. Specifically, if M g is 
reinterpreted as the mass m of a typical star and V s as a 
characteristic relative velocity V* between nearby stars, and 
R 3 is replaced by the typical separation ~ n _1//3 between 
neighbouring stars, one concludes that 



G 2 m 2 n 
V? ' 



(18) 



where n represents a characteristic stellar density. Equation 
(18) agrees (cf. Binney & Tremaine 1987) with the stan- 
dard expression for tn, modulo the absence of the Coulomb 
logarithm In A which would ordinarily appear in the numer- 
ator. This additional factor, which arises in a more careful 
analysis of discreteness effects, reflects the fact that weaker, 
more distant interactions also contribute to energy relax- 
ation, thus reducing the value of ta somewhat. In this sense, 
one infers that an analysis of energy relaxation based on eqs. 
(15) - (17) will, if anything, underestimate the importance 
of environmental effects. 

Estimating the effects of internal pulsations seems more 
uncertain, but one can proceed using dimensional analysis. 
Clearly one can write 



GMg 
tig 



t c ~ PtD, 



and 9 ~ -yV 2 



(19) 



where a, /3, and 7 represent dimensionless constants. Here 
the constant a essentially sets the characteristic amplitude 
of the irregularities associated with the internal pulsations 
as compared with the bulk force associated with the galaxy 
as a whole and, as such, should be small compared with 
unity. One would thus anticipate that a< 1. Appropriate 
values for the constants /3 and 7 are less clear, but it would 
at least seem reasonable to assume that /3 is not signifi- 
cantly larger than unity, i.e., the autocorrelation time t c is 
not much larger than the dynamical time to- The charac- 
teristic size of 7, which enters into D but not n, is more 
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difficult to estimate. Fortunately, however, it follows from 
the simulations described in Section 3 (cf. eq. [12]) that the 
rate of energy relaxation is fixed completely by D and t c , 
independent of r\, so that its precise value is immaterial. In 
any event, given this scaling one infers that 



D 



2„G 2 M 2 q 



2 a 

a p- 



,G 2 M 2 g 



(20) 



where one has used the fact that to ~ R g /V* 



3 NUMERICAL EXPERIMENTS 
3.1 What was computed 

Solutions to the Langevin equation (1) for a constant co- 
efficient of dynamical friction and coloured noise sampling 
the Ornstein-Uhlenbeck process (7) were obtained for three 
different classes of potentials. The most interesting astro- 
nomically, but by far the most expensive computationally, 
were the triaxial generalisations of the Dehnen (1993) po- 
tentials, which have been considered extensively by Merritt 
and collaborators (e.g. Merritt & Fridman 1996). These cor- 
respond to potentials generated self-consistently from the 
triaxial mass density 



p(m) 



with 



(3-7) 
4nabc 



m 7 (1 + m) 



-(4- 7 ) 



x 2 y 2 z 2 
a 2 b 2 c 2 



(21) 



(22) 



for various choices of axis ratios a : b : c and cusp index 
7. Because of the cost associated with integrations of this 
potential, a careful exploration of parameter space was pro- 
hibitive, so that more systematic investigations were per- 
formed instead for the toy potentials 



V(x,y, z) = -(a x +b y + c z ) - 



Mbh 



(23) 



with r 2 — x 2 + y 2 + z 2 , which have been shown (Kandrup & 
Sideris 2000) to reproduce much of the observed behaviour 
for orbit ensembles evolved in the triaxial Dehnen potentials. 
To test the genericity of the basic conclusions, integrations 
were also performed for a completely different class of com- 
putationally inexpensive potentials with 

V(x,y, z) = -(x 2 + y 2 + z 2 ) + ^(x 2 + y 2 + z 2 ) 2 

1 / 2 2, 2 2 . ,2 2x 

— -(x y + ay z + bz x ). 



(24) 



These potentials, which were explored systematically in 
Kandrup, Pogorelov, & Sideris (2000), constitute three- 
dimensional generalisations of the two-dimensional dihedral 
potential first considered by Armbruster, Guckenheimer, & 
Kim (1989). 

The experiments that were performed each involved se- 
lecting ensembles of 8000 initial conditions, all with the 
same energy, and, for given choices of 77, O, and t c , inte- 
grating these initial conditions into the future for a time 
t ~ 100—200 to ■ In order to determine whether energy relax- 
ation impacts regular and chaotic orbits in the same fashion, 
regular and chaotic ensembles were considered separately. 



For the potentials considered, the typical particle speed v 
and the dynamical time to are both of order unity. For this 
reason, most of the experiments assumed that ~ 1. The 
coefficient of dynamical friction rj, which defines the diffu- 
sion constant D = 2O77, was allowed to vary in the range 
10 -8 <f)< 10~ 2 . The autocorrelation time t c was chosen 
to satisfy 10" 2 < t c < 10 3 . It was verified explicitly that in- 
tegrations with autocorrelation times as short as t c = 0.01 
yield results essentially indistinguishable from integrations 
with white noise. 

White noise integrations were performed using a fixed 
time step algorithm developed by Griner, Strittmatter, & 
Honerkamp (1988). Coloured noise integrations were per- 
formed using an extension of this algorithm developed by I. 
V. Pogorelov as part of his Ph. D. thesis (see Pogorelov & 
Kandrup 1999). The basic idea underlying this extension is 
to (i) use a pseudo-random number generator to generate 
white noise X(t), and then (ii) transform X(t) into coloured 
noise Y(t) sampling the Ornstein-Uhlenbeck process by solv- 
ing the stochastic differential equation 



d J + *Y = X(t), 



(25) 



with a — t^ 1 ■ Viewed as an initial value problem (cf. Chan- 
drasekhar 1943b), eq. (25) has the property that any non- 
trivial dependence on initial conditions dies away exponen- 
tially on a time scale a -1 , so that, at late times, the random 
variable Y(t) satisfies 



(Y(t)) = 



and 



<Y(ii)Y(t 2 > = 2CTp(-a|ti-ta|).(26) 



In other words, the late time solution to eq. (25) for X(t) 
given as white noise is coloured noise sampling the Ornstein- 
Uhlenbeck process Y(t). 

3.2 What was found 

Perhaps the single most important conclusion derived from 
these simulations is that the observed behaviour of different 
orbit ensembles is independent of energy and choice of po- 
tential, i.e., the results appear to be universal. Modulo the 
numerical constant A entering into the relation (v 2 ) = AE, 
which is exploited in the derivation of eqs. (9) and (10), 
different ensembles yielded very similar results. This is not 
surprising. When considering the problem of diffusion on 
a constant energy hypersurface, different types of orbits 
and/or different potentials can yield very different behaviour 
(cf. Siopis & Kandrup 2000): Regular orbits are confined 
to three-dimensional phase space tori whereas chaotic or- 
bits move on a higher-dimensional phase space hypersurface 
which is impacted by a complex Arnold web, the details of 
which can depend sensitively on the choice of potential and 
initial condition. However, energy relaxation involves mo- 
tion orthogonal to the hypersurfaces of constant energy, so 
that most of these details are irrelevant. As discussed more 
carefully below, only for very long autocorrelation times t c 
does it appear that the qualitative results depend at all on 
the form of the potential or whether the orbits be regular as 
opposed to chaotic. 

It was also found that, as would be expected, coloured 
noise acts diffusively, so that the root mean squared energy 
for an orbit ensemble grows as the square root of the inte- 
gration time, i.e., SE rms cct 1 ^ 2 . This is, e.g., evident from 
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1 28 
t 



'92 



256 



Figure 1. \SE^ ms \/\Eg\, the fractional change in squared energy 
for an 8000 orbit ensembles evolved for a time t = 256 in the 
generalised dihedral potential with E = 4 and a = b = 1 in the 
presence of coloured noise with © = 1 for several different choices 
of r) and t c . (a) The top curve has r\ = 10 -4 and t c = 10.0. The 
lower has r\ = 10 — 6 and t c = 1.0. (b) The top curve has r\ = 10~ 5 
and t c = 0.1. The lower has 77 = 10 — 5 and t c = 1.0. 
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Figure 3. (a) |<5-E rms |/|i?o|, the fractional change in energy for 
8000 orbit ensembles evolved in the dihedral potential with E = 4 
and a = b = 1 after a time t = 256, viewed as a function of t c for 
different choices of D. The solid line has slope —1.0. 




-8 -6 -4 -2 

logic D 

Figure 2. (a) |(S£ rm3 |/|i?o|, the fractional change in energy for 
8000 orbit ensembles evolved in the dihedral potential with E = 4 
and a = b = 1 after a time t = 256, viewed as a function of D for 
different choices of t c . The solid line has slope 0.5. 

FIG. 1, which exhibits data for an ensemble of 8000 chaotic 
orbits with E = 1 evolved in the generalised dihedral po- 
tential for a total time t = 256 with a — b = 1, allowing for 
coloured noise with = 1 for several different choices of r\ 
and t c . 

Similarly, as for the case of white noise, one finds that 



the individual values of and r\ are irrelevant. For fixed 
autocorrelation time t c , all that matters is the diffusion con- 
stant D — 2Qrj. Moreover, as for the case of white noise, the 
amount of energy relaxation, as probed by 8E rms , scales as 
D 1 ^ 2 for fixed t c . This fact is illustrated in FIG. 2 which, for 
the ensemble of initial conditions used to generate FIG. 1, 
exhibits the fractional root mean squared change in energy, 
|<5-Erm S |/|-Eo|, a t time t = 256 as a function of D for various 
choices of t c . 

For fixed diffusion constant D, the dependence on t c is 
less trivial. As expected, one discovers that, for t c <C ijj, the 
effects of coloured noise are virtually identical to the effects 
of white noise, which can be viewed as a singular limit of 
eq. (5) or (11). Alternatively, for t c ^> tr>, the amount of en- 
ergy relaxation for fixed D is significantly reduced. This is, 
e.g., illustrated in FIG. 3, which, once again for the ensem- 
ble used in FIG. 1, exhibits |5i?rms |/|-Eo| at time t = 256 
as a function of t c for various choices of D. It is evident 
that, even for t c 3> to, the slopes of the curves in this Fig- 
ure are not completely constant, but it does appear that, for 
< l°gio ^ — 2, i.e., for t c between roughly to and lOOtrj, 
the amount of energy relaxation scales as to ■ In this ap- 
proximation, the results illustrated in FIGS. 1-3 yield the 
scaling summarised in eq. (12). 

It is well known that, as far as energy relaxation is con- 
cerned, white noise impacts regular and chaotic orbits iden- 
tically. This also appears to be the case for coloured noise, 
at least for values of the autocorrelation time t c < 30tn or 
so. Interestingly, however, there are distinct differences for 
very large values of t c . As illustrated in FIG. 4, it appears 
that, for t c > 30to or so, regular orbits are impacted more 
than chaotic orbits. In particular, for very large t c the scal- 
ing relation (12) drastically underestimates the importance 
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Figure 4. (a) |<5-E rm3 |/|i?o|, the fractional change in energy for 
8000 orbit ensembles evolved in the dihedral potential with E = 4, 
a = b = 1, and D = 10 — 4 after a time t = 256, for variable 
t c - Triangles refer to an ensemble of regular initial conditions; 
diamonds refer to an ensemble of chaotic initial conditions, (b) 
The same for D = 10 -5 . (c) The same for D = 10~ 6 . 



of energy relaxation. Why this is the case is not completely 
clear. Given, however, that coloured noise acts via a reso- 
nant coupling to the orbits, it seems reasonable to conjecture 
that this effect arises because regular orbits can couple more 
coherently to the slowly varying random forces. 



4 PHYSICAL IMPLICATIONS 

When applied to eq. (16), the expressions for D and t c mo- 
tivated in Section 2.3 allow one to relate the time scale as- 
sociated with energy relaxation to the dynamical time to- 
Consider first the effects of a near-random external en- 



vironment where, presumably, 
infers that 



t c > to- In this case, one 



5E r 



\E\ 



D_ 

\E\ 

G 2 M 2 
B$\E\ 



t 



(27) 



Suppose, however, that the galaxy is at least roughly in virial 
equilibrium. In this case, the 'Virial Theorem' implies that 



GM a /R g ~ \E\ 



(28) 



and the dynamical time to may be taken to satisfy 



i 2 



R'a 



1 



Gp {GM g 



(29) 



Substitution of eqs (28) and (29) into eq. (27) then leads to 
the conclusion that 



5E r 



(R^V (tn\ 1/2 (±V /2 
\RsJ \t c J \t D ) 



(30) 



In a similar fashion, one can allow for the effects of 
internal pulsations, assuming, as in Section 2.3, that t c is 
less than or comparable to to. In this case, 5E rm ,s satisfies 

2 



SE r 



E 



\E\ 



to t, 



(31) 



and by again exploiting eqs. (28) and (29), one concludes 
that 



SE r 



E 



t V' 2 



(32) 



Equations (30) and (32) allow one to infer the typical 
degree of energy relaxation that can be expected over a time 
scale ~ tn, the age of the Universe. 

Consider, e.g., the effects of a dense cluster environment 
where, not unrealistically, the typical velocity dispersion as- 
sociated with stars in an individual galaxy is comparable in 
magnitude to the velocity dispersion of galaxies moving in 
the cluster, and where the typical distance between galaxies 
is of order six times the size of a typical galaxy. In this case, 
R s ~ GR g and t c ~ 6tr>, so that, assuming that tn ~ 100to, 
one would anticipate a fractional root mean squared energy 
8Erms/\E\ ~ 0.1 over the lifetime of the galaxy. Note that, 
as discussed in Section 2.3, this estimate likely understi- 
mates by a small amount the efficacy of energy relaxation 
induced by random interactions with nearby galaxies. 

Similarly, for the case of internal pulsations, it might 
seem plausible to assume that the characteristic fractional 
amplitude of the perturbations a ~ 10 -2 and that t c is com- 
parable to to, so that p ~ 1. One would then anticipate once 
again that, within a time tu ~ lOOto, the fractional root 
mean squared change in energy will have grown to a value 

8Erms/\E\ ~ 0.1. 

These estimates effects - changes in energy by of or- 
der 10% over the age of the Universe - are not huge, but 
they could prove important in real galaxies. Recent numeri- 
cal work suggests that, because of the role of chaos in certain 
key phase space regions, e.g., near corotation and other res- 
onances in spiral galaxies (cf. Wozniak 1993) or in the cen- 
tral regions of cuspy-core ellipticals {cf. Merritt & Fridman 
1996), it may be hard for a galaxy to settle down towards a 
true collisionless equilibrium. However, it would seem sub- 
stantially easier (cf. Siopis & Kandrup 2000) for the system 
to evolve towards a configuration which, albeit not corre- 
sponding to a true equilibrium, could persist in isolation as 
a near-equilibrium for times long compared with the age of 
the Universe. 

The important point, then, is that such near-equilibria, 
even if nearly stable as isolated entities over intervals ~ tn, 
could be destabilised by low amplitude perturbations of the 
form to which real galaxies are necessarily subjected. There 
is compelling numerical evidence that chaotic orbits are 
very susceptible to perturbations which induce phase space 
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transport on (or very nearly on) the constant energy hy- 
persurfaces (cf. Kandrup, Pogorelov, & Sideris 2000, Siopis 
& Kandrup 2000), thus facilitating extrinsic diffusion along 
the Arnold web. Such diffusion will occasion changes in the 
phase space density, which in turn alters the form of the 
gravitational potential, which could in principle triggering a 
nontrivial secular evolution. 

One might, however, argue that this need not have sig- 
nificant implications for bulk evolution since most of the 
phase space is (presumably) dominated by regular orbits 
which (cf. Binney 1978) support the skeleton of the system 
and which, significantly, are not susceptible to such phase 
space diffusion. The key recognition, therefore, is that en- 
ergy relaxation impacts both regular and chaotic orbits; and 
that even comparatively small effects acting on regular or- 
bits, combined with larger effects acting on the (presumably) 
smaller measure of chaotic orbits, could effectively desta- 
bilise a near-equilibrium on a time scale < tu- 

One plausible scenario might involve the effects of in- 
ternal oscillations triggered by a near-collision with another 
galaxy. Numerical simulations (cf. Vesperini & Weinberg 
2000) suggest that such near-collisions can serve to induce 
comparatively long lived, large amplitude oscillations which 
could, e.g., account for the offset nuclei and/or lopsided and 
warped discs observed in certain galaxies (cf. Kornreich, 
Haynes, & Lovelace 1998 and references cited therein); and 
even if these oscillations are dominated by a small number of 
lower order normal modes, they should be accompanied by a 
larger number of higher order modes which could reasonably 
be modeled as near-random excitations. 

In this context, it important to note that energy relax- 
ation induced by near-random internal irregularities is likely 
to be more important for the bulk structure of a galaxy than 
more systematic, nearly periodic effects. Periodic, or near- 
periodic, driving, associated with a small number of com- 
panion objects or with a small number of normal or pseudo- 
normal modes, involves perturbations with a countable set 
of discrete frequencies. The obvious point, then, is that since 
these perturbations have their effect because of a resonant 
coupling to orbits in the galaxy, they can only have a large ef- 
fect on orbits which have considerable power at or near these 
special frequencies and/or harmonics thereof. Chaotic orbits 
typically have broad band Fourier spectra (cf. Tabor 1989) 
and, as such, can couple efficiently to such periodic distur- 
bances (cf. Kandrup, Abernathy, & Bradley 1995). However, 
for regular orbits the power is concentrated at a few special 
frequencies and, if these frequencies are not in a resonance 
or near-resonance with the frequencies of the perturbations, 
the perturbation will have a comparatively minimal effect. 
Periodic driving could perhaps prove important for chaotic 
orbits in general or for regular orbits in specific phase space 
regions. However, one might anticipate generically that it 
will be less important on the overall phase space structure 
which, presumably, is dominated by regular orbits with fre- 
quencies that are not in resonance with the periodic driving. 
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